The symmetry energy at subnuclear densities and nuclei in neutron star crusts 
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We examine how the properties of inhomogeneous nuclear matter at subnuclear densities depend 
on the density dependence of the symmetry energy. Using a macroscopic nuclear model we calculate 
the size and shape of nuclei in neutron star matter at zero temperature in a way dependent on 
the density dependence of the symmetry energy. We find that for smaller symmetry energy at 
subnuclear densities, corresponding to larger density symmetry coefficient L, the charge number of 
nuclei is smaller, and the critical density at which matter with nuclei or bubbles becomes uniform is 
lower. The decrease in the charge number is associated with the dependence of the surface tension 
on the nuclear density and the density of a sea of neutrons, while the decrease in the critical density 



o 

■ can be generally understood in terms of proton clustering instability in uniform matter. 
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Cn ! I- INTRODUCTION 



The inner crust of a neutron star consists of a Coulomb lattice of nuclei embedded in a roughly uniform neutralizing 
background of electrons and in a sea of neutrons 0. The equation of state (EOS) of uniform nuclear matter and the 
surface tension, both at large neutron excess, are essential to understanding of the equilibrium properties of matter in 
the crust. However, laboratory data on nuclei only reflect the bulk and surface properties of nearly symmetric nuclear 
Q^ ' matter So far, calculations of the equilibrium properties of matter in the crust depend on the way to extrapolate 

' the known bulk and surface properties to large neutron excess, which is different among earlier investigations 

In this paper we systematically analyze the question of how the equilibrium properties of inhomogeneous nuclear 
matter at subnuclear densities depend on the parameter characterizing the density dependence of the symmetry energy. 
In doing so, we utilize a macroscopic nuclear model in which the equilibrium nucleon distribution depends on the 
EOS of nuclear matter through minimization of the energy density functional. One of the most important quantities 
is the charge number of the equilibrium nuclide. In a liquid-drop picture 0, this charge number is determined by the 
^ ] size equilibrium condition which controls the ratio between the Coulomb and surface energies. This condition tells 
^ ' that the charge number squared is proportional to the surface tension and the nuclear volume. As we shall see, the 
density dependence of the symmetry energy, which controls the surface tension by affecting the nuclear density and 
the density of the neutron gas, controls the charge number as well. 

We also address the question of how matter with nuclei or bubbles melts into uniform matter with increasing 
density. In this melting process, rodlike and slablike nuclei embedded in a gas of neutrons, as well as rodlike and 
roughly spherical neutron-gas regions (bubbles) surrounded by a nucleon liquid, are expected to occur 
At a density where roughly spherical nuclei are so closely packed that they occupy about 1 /8 of the system volume, 
the nuclei tend to be elongated and eventually fuse into nuclear rods. The advantage of this rod formation is a 
reduction in the total surface area from the roughly spherical case. However, whether bubbles and nonsphcrical nuclei 
actually appear in neutron star crusts depends on the critical density at which proton clustering instability occurs 
in uniform nuclear matter @; they are expected to appear when the density corresponding to the nuclear volume 
fraction of about 1/8 is smaller than the critical density for proton clustering. We find that this critical density is in 
turn controlled by the symmetry energy at subnuclear densities. 

Earlier investigations on such exotic nuclei are more or less based on specific nuclear models P, . An exception 
is the work by Watanabe et al. which is systematic in the sense that the liquid-drop model calculations were 
performed in a way dependent on the proton chemical potential in pure neutron matter, , and the surface tension. It 

was found that the density at which the system dissolves into uniform matter increases with increasing /ip'^''. However, 
it remains to be clarified why some nuclear models [gI IiiIIT^ predict the absence of bubbles and nonsphcrical nuclei. 
It is important to note that these models predict relatively high pressure for pure neutron matter (or, equivalently, 
relatively small symmetry energy) at densities around half the normal nuclear density, while the work by Watanabe 
et al. llQj uses a parametrization jl3l | based on the microscopic calculations by Siemens and Pandharipande as 
the EOS of pure neutron matter and fix the density dependence of the symmetry energy. This parametrization is 
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consistent with the recent Green's function Monte Carlo (GFMC) calculations [T3 at neutron densities up to about 
half the normal nuclear density. We will give a unified picture about the presence of bubbles and nonspherical nuclei 
by describing the pressure of pure neutron matter in terms of the density dependence of the symmetry energy. 

The size and shape of nuclei in the crust bear relevance to the thermal and rotational evolution of neutron stars. This 
is because thermal conductivity and neutrino emissivity in the crust are controlled by electron-nucleus scattering 16j, 
while the motion of superfluid neutron vortices in the crust is affected by vortex-nucleus interactions 17] . However, 
real crustal matter has to be accompanied by defects and impurities, which can play a more important role in the 
star's evolution p^ . In a real neutron star, furthermore, the nuclear system is more or less out of equilibrium in 
the course of mass accretion onto the surface of the star and the star's spin-down [Tgl f20ll2l| . Such disordered and 
nonequilibrium properties are beyond the scope of this paper. 

In Sec. ^ we construct a model for inhomogeneous nuclear matter at subnuclear densities in a way dependent on 
the EOS of nearly symmetric nuclear matter near the saturation point. The equilibrium size and shape of nuclei at 
given density are then calculated from the model constructed in Sec. IIIII Section Hvl is devoted to evaluations of the 
critical density at which uniform matter becomes unstable against proton clustering. Our conclusions are given in 
Sec. El 



II. MODEL FOR MATTER AT SUBNUCLEAR DENSITIES 



In this section, we construct a macroscopic model for zero-temperature, /3-equilibrated, inhomogeneous nuclear 
matter at subnuclear densities. This is an extension of Ref. ^ to the case in which a gas of dripped neutrons is 
present, which is based on Ref. J]. Here we focus on how macroscopic properties of the system depend on the EOS 
of nearly symmetric nuclear matter a nd, for simplicity, ignore various effects such as nucleon pairin g eff ects (2^ , shell 
effects in inhomogeneous matter [23, llj , fluctuation- induced displacements of nuclei and bubbles jlfl , and electron 
screening effects |25l| . 

The bulk energy per nucleon is an essential ingredient of the macroscopic nuclear model. We set this energy as 

3?i2(37r2)2/3 

+ {l ~ a^)vs{n)/n + a^Vn{n)/n, (1) 



where 



and 



am^ + (2) 

1 -I- a^n 



1 b^n^ 

Vn = bin + — — — (3) 
1 + 03/1 

are the potential energy densities for symmetric nuclear matter and pure neutron matter, Un and Up are the neutron 
and proton number densities, n = n„ + rtp, a = {n.n — np)/n is the neutron excess, and to„ is the neutron mass. 
Expressions can well reproduce the microscopic calculations of symmetric nuclear matter and pure neutron 

matter by Friedman and Pandharipande (26|| in the variational method. In this method, the isospin dependence of 
asymmetric matter EOS is shown to be well approximated by Eq. ^ |23| . (Replacement of the proton mass nip by 
rn„ in the proton kinetic energy would make only a negligible difference.) For the later purpose of roughly describing 
the nucleon distribution in a nucleus, we incorporate into the potential energy densities l(2Jl and l(2Jl a low density 
behavior oc n^ as expected from a contact two-nucleon interaction. 

A set of expressions COl-® is one of the simplest that reduces to the usual form Q in the limit of n — > riQ and 



w = u)o + -r^{n - no) + 
18?ig 



So + 77— - "0) 
3no 



(4) 



Here wq, no, and Ko are the saturation energy, saturation density, and incompressibility of symmetric nuclear matter. 
The parameters L and So are associated with the density dependent symmetry energy coefficient S{n): So is the 
symmetry energy coefficient at n = no, and L = 3no{dS / dn)n=no is the symmetry energy density derivative coefficient 
(hereafter referred to as the "density symmetry coefficient"). As the neutron excess increases from zero, the saturation 
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point moves in the density versus energy plane (see, e.g., the dotted hnes in Fig. 2). This movement is determined 
mainly by the parameters L and Sq. Up to second order in a, the saturation energy Ws and density rig are given by 



= Wo + S'oa 



and 



The slope, y, of the saturation line near a = is thus expressed as 



(5) 



(6) 



(7) 



We determine the parameters oi, • • • , 63 in such a way that the charge number, charge radius, and mass of stable 
nuclei calculated in a macroscopic nuclear model constructed in Ref. are consistent with the empirical data. In the 
course of this determination, we fix 63, which controls the EOS of matter for large neutron excess and high density, 
at 1.58632 fm'^. This value was obtained by one of the authors 7} in such a way as to reproduce the neutron matter 
energy of Friedman and Pandharipande 26] . Change in this parameter would make no significant difference in the 
determination of the other parameters and the final phase diagram. 

We describe macroscopic nuclear properties in a way dependent on the EOS parameters ai,---,63 by using a 
Thomas- Fermi model 0. The essential point of this model is to write down the total energy of a nucleus of mass 
number A and charge number Z as a function of the density distributions n„(r) and Tip(r) in the form 



E = Eb + Eg+Ec+ Nninc^ + ZnipC^, 



where 



is the bulk energy. 



Eh 



d^rn{r)w (n„(r), np{r)) 



Eg=Fo I d^r|Vn(r)| 



is the gradient energy with adjustable constant Fq, 

=2 



E^ 



c 



e 



(8) 



(9) 



(10) 



(11) 



is the Coulomb energy, and iV = ^ — Z is the neutron number. This functional allows us to connect the EOS and the 
nuclear binding energy through the bulk energy part Eh- For simplicity we use the following parametrization for the 
nucleon distributions ni{r) (i = n,p): 



ni{r) 



1- ^ 



0, 



r > R, 



(12) 



This parametrization allows for the central density, half-density radius, and surface diffuseness for neutrons and 
protons separately. In order to construct the nuclear model in such a way as to reproduce empirical masses and 
radii of stable nuclei, we first extremized the binding energy with respect to the particle distributions for fixed mass 
number, five EOS parameters, and gradient cocfhcicnt. Next, for various sets of the incompressibility and the density 
symmetry coefficient, we obtained the remaining three EOS parameters and the gradient coefficient by fitting the 
calculated optimal values of charge number, mass excess, root-mean-square (rms) charge radius to empirical data for 
stable nuclei on the smoothed [3 stability line In the range of the parameters < L < 160 MeV and 180 MeV 
< Kq < 360 MeV, as long as ?/ < —200 MeV fm^, we obtained a reasonable fitting to such data (see Fig. 1). As a 
result of this fitting, the parameters no, wq, Sq, and Fq are constrained as no = 0.14-0.17 fm"'^, wq = —16 ± 1 MeV, 
So = 25-40 MeV, and Fq = 66 ± 6 MeV fm^. We remark that a negative L is inconsistent with the fact that the size 
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FIG. 1: (Color online) The sets of (L, Kq) (crosses) consistent with the mass and radius data for stable nuclei. The thin lines 
are lines of constant y. The labels A-I denote the sets for which we perform detailed calculations of the ground state properties 
of inhomogeneous nuclear matter at subnuclear densities. For comparison, the values calculated from two mean-field models 
[TMl (square) and SIII (dot)], which are known to be extreme cases i2Sl . are plotted. The plot shows that our sets of (L, Kq) 
effectively cover such extreme cases. 




oi A = 17, 20, 31 isobars deduced from the experimental values of the interaction cross section tends to increase with 
neutron/proton excess j29l |. This inconsistency can be seen from Eq. ® which shows that the saturation density iig 
increases (and hence the isobar size decreases) with neutron/proton excess for a negative L. 

We remark that in this range the calculations agree well with a more extended data set of nuclear masses for A>2 
IP and charge radii for A > 50 |ai|. The rms deviations of the calculated masses from the measured values are ~ 3-5 
MeV, which are comparable with the deviations obtained from a Weizsacker Bethe formula, while the rms deviations 
of the calculated charge radii from the measured values are about 0.06 fm, which arc comparable with the deviations 
obtained from the A^/'^ law. 

Let us summarize the macroscopic nuclear model used here. This model can describe global nuclear properties such 
as masses and rms radii in a manner that is dependent on the EOS of nuclear matter. One of the important predictions 
of this model was that the matter radii depend appreciably on the density symmetry coefficient L, while being 
almost independent of the incompressibility Kq. Although the present macroscopic approach has some limitations 
in describing the nuclear surface, it is still useful for examining the phase diagram of nuclear matter at subnuclear 
densities 0. 

For the purpose of describing matter in neutron star crusts, we proceed to extend the above-described nuclear 
model to the case of nuclei of various shapes embedded in a gas of dripped neutrons by following a line of arguments 
of Ref. 7]. Here we also take into account a gas of electrons as a constituent of matter in the crust and impose (3 
stability and charge neutrality in the system. 

We consider five phases that consist of spherical nuclei, cylindrical nuclei, planar nuclei, cylindrical bubbles, and 
spherical bubbles, respectively. Each phase is assumed to be composed of a Coulomb lattice of a single species of 
nucleus or bubble at a given baryon density n^. For the convenience of practical calculations, we adopt the Wigner- 
Seitz approximation. In this approximation, a cell in the bcc lattice, including a spherical nucleus or bubble, is 
replaced by a Wigner-Seitz cell defined as a sphere having the same volume (a^) and center. We refer to a as the 
lattice constant. A cylindrical nucleus or bubble having an infinitely long axis and a circular section is taken to be 
contained in a cylindrical Wigner-Seitz cell having the same axis in place of a cell in the two-dimensional triangular 
lattice. For a planar nucleus, a Wigner-Seitz cell is identical with a cell in the one-dimensional layered lattice. For 
the sake of convenience, we redefine the cylindrical and slab Wigner-Seitz cells as a cylinder of height a and base area 
0? and a slab of thickness a and surface area o?, respectively (see Figs. 1 and 2 in Ref. 01). 

For each unit cell, we write the total energy as 

W ^Wn + We + Wc, (13) 

where WM^ We, and Wc are the nuclear energy, the electron energy, and the Coulomb energy. 
As in Eq. (|SJ), the nuclear energy is again expressed in the density functional form: 

HV . / d»,1„(.,»(„„(.,,„,.,r)) 

+mnc^nn{v) + mpC^np{Y) + Fo|Vn(r)n. (14) 
For a spherical nucleus in vacuum, this expression reduces io E — Eq [see Eq. ©]. 
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The electron energy can be approximated as the energy of an ideal uniform Fermi gas, 

^ = ^{^e(2x^ + l){xl + - ln[Xe + {xl + 1)^']} (15) 

with 

Xe ; 

rrieC 

where is the electron mass, and rig is the electron number density that satisfies the charge neutrality condition, 



(16) 



a^He = / (frnp{r). (17) 

We remark that rtg is so high that we can safely ignore inhomogeneity of the electron density induced by the electron 
screening of nuclei or bubbles and the Hartree-Fock corrections to the electron energy. 

The Coulomb energy is composed of the proton self-Coulomb energy and the lattice energy. We write the Coulomb 
energy as 



Wc^l- I d're[np{r) - nM^^) + AM^i, (18) 

^ ./noli 



where 0(r) is the electrostatic potential in a Wigner-Seitz cell, and AWi is the difference of the rigorous calculation 
[3^ for a cell in the bcc (triangular) lattice of spherical (cylindrical) nuclei or bubbles having sharp surfaces from 
the Wigner-Seitz value, as parametrized in Ref. "T"!. We take into account AWi, which is a less than 1 % correction, 
because AWi depends sensitively on the dimensionality of the lattice. (Note that AWi = for the layered lattice of 
slab nuclei.) 

For nucleon distributions in the Wigner-Seitz cell, we simply generalize the parametrization (|12() for a nucleus in 
vacuum into 



, ti 

r 



Ri 



r < Ri, 

(19) 

R^ < r. 



Here r is the distance from the central point, axis, or plane of the unit cell. In the case of nuclei, = 0, while in 
the case of bubbles, n^p —Q. 

We finally determine the equilibrium configuration of the system at given baryon density, 

Ub = a'^ / dVn(r). (20) 

Jccll 

First, for each of the five inhomogcneous phases, we minimize the total energy density W/a^ with respect to the 
eight parameters a, n™, jt-™*, n^p (for nuclei) or (for bubbles), R„, Rp, tn, and tp. This minimization implicitly 
allows for the stability of the nuclear matter region (the region containing protons) with respect to change in the 
size, neutron drip, /3-decay, and pressurization. In addition to the five inhomogeneous phases, we consider a uniform 
phase of /3-equilibrated, neutral nuclear matter. The energy density of this phase is the sum of the nucleon part 
nw + m„c^n„ -|- rUpC^np [see Eq. (QJ] and the electron part l|15fl . By comparing the resultant six energy densities, we 
can determine the equilibrium phase. 



III. EQUILIBRIUM SIZE AND SHAPE OF NUCLEI 

We proceed to show the results for the equilibrium nuclear matter configuration obtained for various sets of the 
EOS parameters L and Kq as shown in Fig. 1. These parameters are still uncertain since they are little constrained 
from the mass and radius data for stable nuclei 01 ■ As we shall see, the charge number of spherical nuclei and the 
density region containing bubbles and nonspherical nuclei have a strong correlation with L. 

We first focus on spherical nuclei, which constitute an equilibrium state in the low density region. We calculate the 
charge number of the equilibrium nuclide as a function of rn, for the EOS models A-I as depicted in Fig. 2. Note that 
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FIG. 2: Energy per nucleon of nuclear matter for the nine sets of {L,Ko) referred to as A-I in Fig. 1. In each panel, the solid 
lines are the energy at a = 0, 0.3, 0.5, 0.8, 1, and the dotted line is the saturation line. 
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FIG. 3: (Color online) The charge number of spherical nuclei as a function of n^, calculated for the EOS models A-I. 



the recent GFMC calculations of the energy of neutron matter based on the Argonne v8' potential [15] are close to 
the behavior of the model E. Hereafter we will thus call the model E as a typical one. The result is shown in Fig. 3. 
For densities below ^ 0.01 fm~^, the calculated density dependence of the charge number Z is almost flat, a feature 
consistent with the results in earlier investigations 0. More important, the calculated charge number is larger for 
the EOS models having smaller L, and this difference in Z is more remarkable at higher densities. 

As we shall see later in Fig. 6, this property of Z is related to the tendency that with increasing L, the nuclear density 
decreases while the density of the neutron gas increases. Note that Z is, within a liquid-drop model jljj, determined by 
the size equilibrium condition relating the Coulomb and surface energies in such a way that Z increases with increasing 
surface tension. Since the Thomas-Fermi model adopted here can be mapped onto a compressible liquid-drop model 
PI, the present results may well be interpreted in terms of the liquid-drop model. In fact we shall estimate the surface 
tension from the Thomas-Fermi model as a function of L and discuss how the surface tension depends on the nuclear 
density and the neutron sea density. 

We also note that the density at which the phase with spherical nuclei ceases to be in the ground state is between 
0.05 fm~'^ and 0.07 fm^"^. This result, consistent with the results obtained in earlier investigations 0,0,^3' '^il^ be 
discussed below in terms of fission instability. 

The average proton fraction, which is the charge number divided by the total nucleon number in the cell, is plotted 
in Fig. 4. We observe that the dependence of the average proton fraction on the EOS models is similar to that of Z. 
We also find that the average proton fraction basically decreases with baryon density. This is a feature coming from 
the fact that as the baryon density increases, the electron chemical potential increases under charge neutrality and 



7 




Hb (fm 

FIG. 4: (Color online) The average proton fraction as a function of rib, calculated for the EOS models A-I. 
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FIG. 5: (Color online) The density region containing bubbles and nonspherical nuclei as a function of L, calculated for the EOS 
models A-I. For comparison, the density corresponding to it = 1/8 in the phase with spherical nuclei and the onset density, 
n{Q), of proton clustering in uniform nuclear matter, which will be discussed in Sec. II VI are also plotted by circles and crosses, 
respectively. 



then the nuclei become more neutron-rich under weak equilibrium. 

We next consider the density region where bubbles and nonspherical nuclei appear in equilibrium, i.e., the density 
region of the "pasta" phases. We start with such a density region calculated for the EOS models A-I. The results 
are plotted in Fig. 5. Except for the model C, we obtain the successive first order transitions with increasing density: 
sphere — > cylinder — > slab — > cylindrical hole —>■ spherical hole — > uniform matter. A marked correlation of the upper 
end of the density region with the parameter L can be observed by referring to Fig. 1. This dependence will be 
examined in detail in the next section. 

As can be seen in Fig. 5, the lower end of the density region of the pasta phases has only a weak dependence on the 
EOS models. In order to have a closer look at this feature, it is instructive to calculate a density at which spherical 
nuclei become susceptible to fission-inducing quadrupolar deformations. Within the framework of a liquid-drop model, 
a spherical liquid-drop in a Wigner-Seitz cell is predicted to undergo such a fission-like instability when the volume 
fraction u, i.e., the ratio of the liquid-drop volume to the cell volume, becomes approximately 1/8 1]. In such a 
closely packed situation, the Coulomb self-energy of the liquid-drop amounts to twice the surface energy even under 
size equilibrium. We note that the density corresponding to u = 1/8 is generally within ±0.01 fm^'^ of the transition 
density calculated from the energy comparison between the phases with spherical nuclei and with cylindrical nuclei. 

In the present model, we evaluate the volume fraction u as 47r(rp/a)^/3, where rp is the rms radius of the proton 
distribution multiplied by a factor -y/S/S. This is because the proton self-Coulomb energy is relevant to the fission-like 
instability. At u = 1/8, we calculate the equilibrium properties of matter with spherical nuclei for the parameter sets 
{L,Ko) included in Fig. 1. The results are plotted as a function of L in Fig. 6. The results for L > 100 MeV are 
scarce since in this case the pressure of neutron matter is too high for u to amount to 1/8. We remark that the results 
show only a weak dependence on Kq. 

It is important to note that the baryon density at w = 1/8 is almost flat at ^ 0.06 fm^"^ [see Fig. 6(a)]. This is 
consistent with the lower end of the density region of the pasta phases as depicted in Fig. 5. This magnitude of nt 
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FIG. 6: (Color online) The equilibrium properties of matter with spherical nuclei as a function of L, calculated at fixed volume 
fraction it — 1/8. The baryon density rij, (a), the nucleon densities in the center and boundary of the cell (b), the nuclear 
charge number Z (c), the central and average proton fractions (d), the lattice constant a (e), and the effective surface tension 
a (f) are plotted. 



at u = 1/8 can be roughly understood from a simple formula Ub — (n™ + n™)u + ?t.°"*(1 — u) with the values of njf , 
n°^*, and n°"* in Fig. 6(b). We also note that with increasing L, the central density decreases while the neutron sea 
density increases. This is natural because both the saturation density of nuclear matter at nonzero neutron excess 
and the symmetry energy at subnuclear densities decreases with L. 

We now turn to the L dependence of the charge number at u = 1/8 [see Fig. 6(c)]. The charge number decreases 
with L. This feature can be understood from the size equilibrium condition within a liquid-drop picture 0. This 
condition states that the Coulomb energy of a cell is half as large as the nuclear surface energy. Consequently, the 
equilibrium charge number squared is proportional to the surface tension and to the nuclear volume. Note that 
the proton fraction in the nuclear center and the nuclear volume have a relatively weak dependence on L [see Figs. 
6(d) and 6(e)]. The surface tension is thus expected to have a similar L dependence to that of the charge number 
squared, through the dependence on the densities inside and outside the nucleus. From the Thomas- Fermi model, it 
is reasonable to estimate the effective surface tension as 



(21) 
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where 

Wg= f dVFo|Vn(r)p, (22) 

Jccll 

is the gradient energy per ceU. This is because in equihbrium, the Coulomb energy of a ceU is as large as Wg 1^, 
implying that the nuclear surface energy, w Anarp, is twice as large as Wg. The surface tension thus estimated 
basically follows the behavior of Z^, as can be seen from Figs. 6(c) and 6(f). The surface tension is generally the 
function of the neutron excess in the nuclear interior and the densities inside and outside the nucleus Since the 
density gradient in the surface region tends to become small for smaller difference between the central density and 
the neutron sea density, the surface tension decreases with L as shown in Fig. 6(f). 

IV. PROTON CLUSTERING IN UNIFORM MATTER 

In this section, we focus on the upper end of the density region of the pasta phases. This upper end corresponds 
roughly to a density at which uniform nuclear matter neutralized and /3 equilibrated by electrons becomes unstable 
against proton clustering. In fact, this correspondence can be seen from Fig. 5. 

We calculate the onset density of proton clustering by following a line of argument of Baym, Bethe, and Pethick 
|l3j| . This density was obtained in Ref. TS] by expanding the energy density functional E[ni{r)] {i = n,p,e) of the 
system with respect to small density fluctuations 6ni{r) around the homogeneous state. While the contribution of first 
order in dni(r) vanishes due to equilibrium of the unperturbed homogeneous system, the second order contribution 
can be described in the spirit of the Thomas-Fermi model used here as 

E-E, = \j -^v{q)\5n,{ci)\\ (23) 

where Eq is the ground-state energy, Jrip(q) is the Fourier transform of (5np(r), and v{q) is the potential of the effective 
interaction between protons as given by 



47re 

V{q) = v, + I3q^+ . (24) 



Here, 



dup dun/dun 
P^2Fo{l + 2C + e), (26) 



(27) 



1 /3 

with fj,n{p) the neutron (proton) chemical potential, and fcxr ~ O.Srie is the inverse of the Thomas-Fermi screening 
length of the electron gas. The effective potential v{q) takes a minimum value Vmin at q = Q, where 



(28) 



^^min = ^'o + 2(4^e2/3)V2 _ pklp. (29) 

In the energy expansion up to second order in Srii , the condition that uniform nuclear matter becomes unstable with 
respect to proton clustering reads Vmin = 0. Generally, Wmin is dominated by the bulk contribution vg, which decreases 
with decreasing density (see Fig. 7). This density dependence ensures the presence of a critical density, n(Q), above 
(below) which the matter is stable (unstable) with respect to proton clustering. Hereafter we will estimate n{Q) 
without including the gradient and Coulomb contributions to Wminj which act to reduce n{Q) only by an amount of 
order 0.02 fm^^ (ggg pjg^ 7 ^Iso Ref. 33]). 
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The density dependence of vq can be seen by substituting Eq. into Eq. p5() . In the Umit of a — > 1, to which 
nuclear matter (3 equihbrated and neutrahzed by the electron gas is close at subnuclear densities (see Fig. 7), vq 
behaves roughly as 

^ {d^i/dn)a^o S{n) 

^0 ~ a , (60) 

[ofi/an)a=i n 

with the symmetry energy coefficient S{n) = Wa=i(n) — Wa=o{n-)- At subnuclear densities, S{n)/n depends 
only weakly on n, while the compressibility ratio between pure neutron matter and symmetric nuclear matter, 
{dfi/dn)a=o/{diJ,/dn)a=i, increases almost linearly with density because of the saturation property of the symmetric 
nuclear matter. 

In Fig. 8 we show the results for n{Q) estimated for the parameter sets {L,Ko) included in Fig. 1. We find that 
n{Q) decreases with L, whereas it does not have a marked dependence on Kq. The L dependence is correlated with 
the L dependence of the symmetry energy coefficient S{n) since S(n) acts as a driving force of proton clustering. 
Note the general tendency that at subnuclear densities, the larger L, the smaller symmetry energy coefficient S{n) 
(see Fig. 7). The proton clustering thus takes place at lower density for larger L. Figure 7 also shows that for n > 0.1 
fm~'^, S{n) becomes larger for larger L. This is a feature coming from the empirical relation, Sq « 0.075L + 28 MeV, 
derived in Ref. 

We can also observe from Fig. 8 that the difference between n{Q) and the density corresponding to u = 1/8 in 
the phase with spherical nuclei decreases with L and eventually vanishes near L = 100 MeV. This suggests that the 
density regime of the pasta phases is limited for a large L, although for the standard EOS model E, corresponding 
to L ~ 40 MeV, it does appear between ~0.06 and ^^0.09 fm~^. In our EOS model, a larger value of L implies a 
harder EOS of pure neutron matter as we shall see below. We thus conclude that the absence of the pasta phases 
seen in Refs. |^ from the EOS model with relatively high pressure of neutron matter at subnuclear densities 

is consistent with our result. 

In order to clarify this consistency, we calculate the pressure of pure neutron matter, 



2 

OUn 



(31) 
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FIG. 8: (Color online) The onset density of proton clustering in uniform nuclear matter as a function of L. For comparison, 
we plot the density corresponding to = 1/8 in the phase with spherical nuclei, which is taken from Fig. 6(a). 
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FIG. 9: The proton chemical potential and pressure in pure neutron matter of density 0.1 fm ^ as a function of L. 

for the parameter sets [L^Kq) shown in Fig. 1. The results for P„ at n„ = 0.1 fm~'^ are plotted as a function of 
L in Fig. 9. We find out a roughly linear L dependence of Pn at n„ = 0.1 fm~'^. This dependence can be roughly 
understood by substituting the expansion Q into Eq. 131f) and thereby obtaining 



~?r — ~ "oj + ^"-0 — 

% \na J 3 V"o, 



(32) 



We remark that this pressure controls the neutron skin thickness of ■^"^Pb evaluated within the framework of the 
Skyrme Hartree-Fock model |3^ . 

We conclude this section by mentioning a relation between the systematic liquid-drop analysis [lo| and the present 
analysis. In Ref. the value of L was fixed at L = 60 MeV, and the proton chemical potential in pure neutron 
matter, 



4°^ 



dijiw) 



drir, 



(33) 



was changed by a factor of 2. In the present analysis, on the other hand, the value of L was taken between and 160 
MeV, while /ip"^ depends only weakly on the value of L, as shown in Fig. 9 in which the results for /Zp"' calculated at 
n„ = 0.1 fm~^ for the parameter sets {L,Kq) shown in Fig. 1 are plotted as a function of L. According to Ref. [lol |. 

fi^^ plays a role in shifting the density region of the pasta phases without changing its width significantly. We may 
thus conclude that it is the parameter L that controls the presence of the pasta phases. 



V. CONCLUSIONS 



We have analyzed the equilibrium properties of inhomogeneous nuclear matter at subnuclear densities in a way 
dependent on the density symmetry coefficient L by using a macroscopic nuclear model. We have estimated the upper 
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and lower ends of the density region of the pasta phases from the onset densities of proton clustering in uniform nuclear 
matter and fission-like instability of spherical nuclei, respectively. We find that the upper end decreases with L, while 
the lower end is almost flat at 0.05-0.07 fm^^. The former arises from the L dependence of the symmetry energy, 
while the latter can be understood from the volume fraction u ~ 1/8 at which spherical nuclei become susceptible 
to fission-inducing deformations. For a typical EOS model consistent with the GFMC calculations of pure neutron 
matter, the calculated pasta regime is appreciable. In fact, the pasta regime is predicted to appear when L ^ 100 
MeV. 

The present analysis is the first to attempt a systematic analysis of the pasta region in terms of L. However, much 
care needs to be taken of the interpretation of the results. While L is the parameter characterizing the expansion of 
w with respect to n and a around n = Uq and a = 0, the system of interest here is at large neutron excess and at 
subnuclear densities. The relation between the parameter L and neutron star matter depends on the parametrization 
of w with respect to n and a. It is thus useful to keep in mind that we confined ourselves to expressions C3)^® 
although they are known to be capable of reproducing various existing microscopic calculations of the EOS of uniform 
nuclear matter. 

We have also calculated the charge number Z of spherical nuclei as a function of density for various values of L. 
Generally, the charge number Z becomes smaller for larger L, a feature that could be of relevance to the evolution 
of neutron stars |l6j|. In order to make better estimate of Z, however, shell and pairing effects should be taken into 
account. 
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